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ABSTRACT 

Aims. We study the circumstellar evolution of the binary HD 101584, consisting of a post-AGB star and a low-mass companion, which 
is most likely a post-common-envelope-evolution system. 

Methods. We used ALMA observations of the ^^CO, ^^CO, and C^^O / = 2-l lines and the 1.3 mm continuum to determine the 
morphology, kinematics, masses, and energetics of the circumstellar environment. 

Results. The circumstellar medium has a bipolar hour-glass structure, seen almost pole-on, formed by an energetic jet, ^ 150 km s“^ 
We conjecture that the circumstellar morphology is related to an event that took place ^ 500 yr ago, possibly a capture event where the 
companion spiraled in towards the AGB star. However, the kinetic energy of the accelerated gas exceeds the released orbital energy, 
and, taking into account the expected energy transfer efficiency of the process, the observed phenomenon does not match current 
common-envelope scenarios. This suggests that another process must augment, or even dominate, the ejection process. A significant 
amount of material resides in an unresolved region, presumably in the equatorial plane of the binary system. 
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1. Introduction 


The source HD 101584 (V885 Cen, IRAS 11385-5517) is bright 
at optical wavelengths (V ^ 7 mag ) with a spectral type A6Ia 
(ISivarani et all 1 199^ lKippei]l200^ . iBakker (Il996bh pre¬ 
sented optical and infrared data that indicate that it evolved from 
the asymptotic giant branch (AGB) at most a few 100 years ago, 
and they estimated a (present-day) mass of 0.6 M©, a luminos¬ 
ity of 4000 Lq, and a distance of 0.7 kpc. 

Photometric and radial veloci ty variations show that 
HD 10158 4 has a binary companion. iBakker et~^ (Il996ah and 
Diaz et al.] (l2007h f o und pe riods of 218 and 144 days, respec- 
tivelv. IBakker et al.l (Il996ah inferred an essentially edge-on cir- 
cumbinary disk to explain optical absorption lines. The absence 
of spectroscopic emission from the companion indicates that this 
is a low-mass main-sequence star. 

Images from the Hubble Space Telescope (HST) show a dif¬ 
fuse circumstellar environment. T here may be a rin g of radius 
1'.'5 roughly centred on the star (ISahai et al.ll2007l) . that is al¬ 
most circular, which suggests that we see it almost pole-on. The 
central star is bright, presumably because material around the 
polar axis of the system has partially been evacuated. 

The c ircumstellar gas has been detected in molecular line 
emission. lOlofsson & NvmanI (Il999h obtained and 
/ = 1-0 and 2-1 single-dish map data. They found at least three 
kinematic components: a low-velocity, an intermediate-velocity, 
and a high-velocity component of total widths 15, 100, and 

j^300kms“\ respectively. The former two were spatially un¬ 
resolved, and the high-velocity gas shows an east-west bipolar 


outflow with a Hubble-like velocity gradient, the most blue- and 
red-shifted emission at 5" to the W and the E. 

A narrower (^80kms ~^), double-peaked OH 1667 MHz 
maser line was mapped bv Ite Lintel Hekkert et al.l (Il992h . The 
integrated OH emission is centred on the star (within 073), but 
is absent inside a radius of 1". The velocity increases system¬ 
atically along the position angle (PA) -60° with the most blue- 
and red-shifted emission at ^ 7!' to the SE and the NW. OH emis¬ 
sion and a 10 pm feature (IBakker et al.lll996ah indicate an 0-rich 
(C/0<1) circumstellar medium. 

Hence, the circumstellar environment of HD 101584 appears 
to be consistent with that of a classic proto-planetary nebula 
(proto-PN): an equatorial density enhancei nent with an orthog¬ 
onal, bipolar, high-velocity outflow (^e.g.. lAlcoleaet^ 1200 ll 


120071: l^stro-Carrizo et al.ll20^ISahai et al.ll200'^ 

In this Letter we present data acquired with the Atacama 
Large Millimeter/sub-millimeter Arra}Q (ALMA) for the ^^CO, 
^^CO, and C^^O / = 2-l line and the 1.3 mm continuum on 
HD 101584. In addition to this, we have detections of line emis¬ 
sion from SiO, CS, H 2 S, SO, SO 2 , OCS, H 2 CO, and some of 
their isotopologues, which will be presented in forthcoming pub¬ 
lications. 


^ This Letter makes use of the following ALMA data: 
ADS/JAO.ALMA#2012.1.00248.S. ALMA is a partnership of 
ESO (representing its member states), NSF (USA) and NINS (Japan), 
together with NRC (Canada) and NSC and ASIAA (Taiwan), in 
cooperation with the Republic of Chile. The Joint ALMA Observatory 
is operated by ESO, AUI/NRAO and NAOJ. 
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Table 1. Observing parameters 


Tuning 

y 

[GHz] 

^on-source 

[min] 

AV 

[kms“^] 

channel rms 
[mJy beam“^] 

Beam 

["X",o ] 

Flux cal. 

J1107-4449 
Flux [mJy] 

J1131-5818 
Flux [mJy] 


230.53 

8 

0.64 

1.8 

0.65x0.57,8.2 

Ceres 

710 + 0.7 

92.4 + 0.3 

‘^co, 

220.06 

26 

0.67 

0.8 

0.65x0.54,46.1 

Titan 

680 + 2 

86.3 + 0.3 


2. Observations 

The ^^CO, ^^CO, and J = 2-l line data were obtained on 
29 April 2014 with 35 antennas of the ALMA 12 m array in two 
frequency settings in band 6 (cycle 1). In both cases, the data 
contain four spectral windows (spw) of 1.875 GHz width and 
with 3840 channels. The observing parameters of the spw rele¬ 
vant for the data presented here are listed in Table [T] 

The baselines range from 16 to 497 m. This means a maxi¬ 
mum recoverable scale of 10". Bandpass and gain calibration 
were performed on J1107-4449 and J1131-5818 (3° separation 
from the target). Flux calibration was made using Ceres and Ti¬ 
tan for the different tunings (Butler-JPL-Horizons 2012, ALMA 
Memo 594). Based on the calibrator fluxes determined in the two 
tunings, we estimate the absolute flux calibration to be accurate 
to within 5%. 

The data were reduced using CASA 4.2.2. After corrections 
for the time and frequency dependence of the system temper¬ 
atures and rapid atmospheric variations at each antenna using 
WVR data, bandpass and gain calibration were made. Poor phase 
coherence led to the flagging of three antennas in the tun¬ 
ing and one antenna in the tuning. Subsequently, for each 
individual tuning, self-calibration was performed on the strong 
continuum. Imaging was made using the CASA clean algo¬ 
rithm after a continuum subtraction was performed on the emis¬ 
sion line data. Our final line images were created using natural 
weighting and the restoring beam given in Table [T] We produced 
a single continuum map combining both tunings with a line-free 
equivalent bandwidth of 10.2 GHz. 

3. Observational results and discussion 

3.1. Global characteristics 

The 1.3 mm continuum map, overlaid on the global ^^CO(2-l) 
map, is shown in Fig.[T] The total flux density, S 230 = 112mJy, 
falls short by almost a factor of four compared to the extrapo¬ 
lated sp ectral energy distribution (SED) model of iBakker et al.l 
(Il996bl) . suggesting a significant loss of flux originating from a 
region larger than 10". The morphology is complex with faint 
structures surrounding an intense inner region, where 70% of 
the total flux comes from inside a central 2" diameter circle. 

The ALMA global CO line profiles are shown in Fig.|^ The 
systemic velocity, i^sys, is estimated to be 42+ lkms“^ (LSR). 
The ALMA total fluxes, integrated over the global line pro¬ 
files, are 690, 150, and 18Jykms“^ for the ^^CO, ^^CO, and 
/ = 2-l lines, respectively. The isotopologue line inten¬ 
sity ratios in the full velocity range and the central 20 km s“^ are 
12C0/1^C0^5 and 2 (solar value 89) and and 

4 (solar value 6), respectively. This indicates that opacity effects 
are significant in the line, affect the line, and that they 
increase significantly towards the central velocity range. Assum¬ 
ing that the ratio is solar (reasonable for a lower-mass 

star), the may not be too far from solar. 

The single-dish ^^C O(2-l) flux is 850Jykms“^ 
dOlofsson & Nvma^ Il999l) . with an uncertainty of ^^20%, 
that is, some flux is probably missed by ALMA. Most likely. 


Fig. 1. 1.3 mm continuum image (black contours) of HD 101584 su¬ 
perposed on the global ^^CO(2-l) map (synthesized beam in the lower 
left corner). The dynamic-range-limited rms is 0.14mJy beam“^ in a 
beam of 0'.'56x0'.'49 (PA = 41.4°) in the continuum. The contour levels 
are (2, 4, 8, 20, 40, 60, 80, 100) xrms. The line intensity scale unit is 
Jykms“^ beam“k 



Vlsr [km/s] 

Fig. 2. SEST (black, offset +5Jy) and the ALMA global 
(blue, offset +2Jy), (red, offset +lJy, scaled xl.5), and 
(green, scaled x3) / = 2-1 lines. 


based on a detailed comparison of the ALMA global and 
the single-dish line profiles, no flux is missed in the central 
20 km s“^ and at the extreme velocity peaks, while 25% of the 
flux is missed in the velocity range 10 < | v - v^ys I < 90kms“\ 
Fig.|2l This suggests the presence of extended gas that has been 
efficiently accelerated. 

The channel m aps confirm the conclusion by 

lOlofsson & NvmanI (1 19991) of an E-W oriented high-velocity 
outflow, but also show additional considerable complexity. 
Fig. [ 3 ] The morphology of the emission is very symmetric with 
respect to the systemic velocity. The high-velocity outflow is nar¬ 
row and directed at PA 90°. A complex ring-like structure is 
present in the velocity range -30 -120 km s“^. The strong central 
emission is confined to a narrow velocity range, < 20 km s“^ 



4 2 0 -2 -4 

R.A. offset [arcsec] 
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Fig. 3. Velocity-channel maps (Az; = 12.5 km s"^) of the ^^CO(/ = 2-1) 
emission towards HD 101584 (synthesized beam in the lower left cor¬ 
ner). The intensity scale unit is Jybeam"^ (noise rms 2mJybeam"^). 



-100 0 100 200 


LSR velocity [km/s] 

Fig. 4. ^^CO(/ = 2-l) position-velocity diagram along PA = 90°. The 
systemic velocity (solid) and the velocity ranges discussed in text 
(hashed) are indicated. The intensity scale unit is Jybeam“^ 


3.2. A jet, an hour-glass structure, and the central region 

A number of different structures can be identified in different 
velocity ranges. 

\v-Vsy^\> 60kms~^: The position-velocity (PV) di¬ 

agram along the PA = 90° shows a jet-like structure, with 
a Hubble-like velocity gradient, reaching a peak velocity of 
J^150kms“^ Fig. IH The emission is patchy along the jet, 
which ends in two sharp features at the extreme velocities. The 
PNe BD-r30°3639 pr ovides another example of such features 
(iBachiller et al.l[200Ql) . These may indicate that this is as far as 
the jet has travelled since its formation, or it has been stopped 
by surrounding material, or, at these distances, the jet is directed 
along the line of sight as a result of precessing, for example. 

I u-Vsys I <30kms~^: The CO emission in the velocity range 
L'sys ± 30kms“^ is dominantly ring-like, as exemplified by the 
data in Fig. [5] The ring size increases with offset from the 
systemic velocity, and its centre position depends on velocity in 
a manner consistent with the direction of the velocity gradient of 
the high-velocity outflow. The PV diagram. Fig. IH shows 
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Fig. 5. Velocity-channel maps (Au = 4kms"^) of the ^^CO(/ = 2-l) 
emission towards HD 101584 (synthesized beam in the lower left cor¬ 
ner). The intensity scale unit is Jy beam"^ (noise rms 2mJy beam"^). 


the expected morphology of a hollow structure, slightly tilted 
w.r.t. the line of sight at PA 90° in this velocity range, sug¬ 
gesting an hour-glass structure seen almost pole-on. We estimate 
the inclination angle, /, to be about half of the opening angle of 
the hour-glass structure, 10°. Surprisingly, the ring structure is 
elongated along PA 90°, as opposed to the expected PA 0°. 

At lower brightness level there is evidence for another bipo¬ 
lar structure in this velocity range. It has a velocity gradient op¬ 
posite to that of the major CO outflow and a different position 
angle (PAj^- 45°), Figs. [3] and IH The direction of the structure 
and its velocity gradient are roughly consistent with those of the 
OH 1667 MHz line. 

30<\v-Vsys\<60kms~^\ At j^+30kms“^ from i^sys, the 
ring-like emission breaks up into a complex structure, yet very 
symmetric w.r.t. the centre. Fig. O The structure ends at offset 
velocities of j^60kms“^ in two compact features displaced by 
^1" on either side of the centre, along PA 90°. Assuming that 
this structure has its origin in a region close to the binary system, 
and adopting a tilt angle of 10°, a projected expansion length of 
1'.'5, and an expansion velocity of 60 km s“^ the age is estimated 
to be 500 yr. With the same assumption on the tilt angle, the 
time scale of the jet is 520 yr. Thus, the complex structure and 
the jet could have their origins in the same event. 

I u-v^ys \< 10kms~^\ Using the 0(2-1) line as a mass in¬ 
dicator, it is clear that a substantial fraction of the mass lies in 
the central region. About 55% of the integrated flux comes from 
the range L'sys ± 10 km s“^. Further evidence for substantial mass 
in this region is the strong and narrow H 2 S( 22 o- 2 ii) line, and the 
corresponding lines from its and isotopologues. Part of 
the emission in this velocity range comes from a circular region 
2" in diameter (possibly related to the ring in the HST im¬ 
ages), and part of the emission comes from an unresolved com¬ 
ponent at the center. The unresolved component is characterized 
by narrow spectral features, only a few km s“^ wide. It is promi¬ 
nent in high-energy 100 K) SO 2 lines, for instance. 


3.3. Mass, energy, and momentum estimates 


A simple estimate of the dust mass using 


Md 


S 23oD^ 
/<23oB230 


( 1 ) 


where B 230 is the blackbody brightness at 230 GHz, with 
^230 = 112mJy, D = 0.7 kpc, the dust opacity /C 230 = 0.5 cm^ g“^ 
(iLiseau etaH 12015), and an adopted du st temperature of 45 K 
[from the SED in Bakker et"^ (Il996bh l gives Md ~ 0.007 Mq 
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(with an uncertainty of, at least, a factor of four). This suggests 
about a solar mass of gas in the central region, assuming a dust- 
to-gas ratio of 200. There probably is cons iderably more mass in 
an extended region iBakker et al.l (Il996bh estimated 0.02 M© in 
total]. 

A simple estimate of the gas mass is obtained using 




hcguA^if 




( 2 ) 


where normal symbols are used for the constants, and / is 
the abundance with respect to H 2 , I the flux density inte¬ 
grated over velocity range, Q the partition function, T^ot the 
rotational temperature, and the upper level energy. Us¬ 
ing the C^^O(2-l) line, /=18Jykms“^ /ci8o=4xl0“^ [so- 
1 i 6 q/ 18 q^ fco^ 2 x l 0 -^ ( 0 -rich envelope)], and T,ot=30K 

(iBujarrabal et al .1120011) . we obtain Mg ^03 Mq (with some con¬ 
siderable uncertainty, e.g., the gas temperature probably varies 
over the source). More than half of it lies in the central 20 km s“^ 
range. 

Using the 0(2-1) line and the same assumptions as 
above, we And that the kinetic energy in the material accelerated 
to velocities > 10 km s"Ms 4 x 10"^^ erg (a lower limit since ex¬ 
tended flux from accelerate d gas is missed). This li es in the mid¬ 
dle of the range found bv IBujarrabal et al.l (l 200 ll) . erg, 

for the bipolar outflows of a sample of 28 proto-PNe. 

Likewise, we And a momentum of j^lO^^gcms"^ in 
the accele rated gas, which lies i n the middle of the range 
found by IBujarrabal et al.l (l 200 ll) . 10 ^^“"^^ gems The es¬ 
timated momentum rate over 500 yr is 7 x 10^^ ergcm“\ 
which is substantially higher than that available from radiation 

I L/c), ^ 5 X 10^^ erg cm ~\ a common situation for proto-PNe 
Bujarrabal et al.|[200lh . Hence, a different momentum source is 
required. 


3.4. A scenario 


Adopting the orbital parameters (period P = 218 days and ve¬ 
locity semi-amplitude ^i=3kms“^) and the mass estima te. 
Ml = 0.6 Mq, of the post-AGB star (iBakker et al.l ll996allEl . 
/ = 10 °, and a circular orbit (reasonable if we have a capture 
event), we derive a mass function of 


(M 2 sin i)^ 
(Ml -r M2)2 


47i^(ai sin i)^ 
GP2 


= 6x lO-^Mo 


( 3 ) 


for the binary system, where M 2 is the mass of the companion 
(ai = KiP/lnsin i). Thus, we estimate that M 2 ~ 0.6 M© at a dis¬ 
tance of a^O.7 A\J (a 2 = aiM 2 lMi, a = ai-\-a 2 ) to the primary. 
This could be the result of a capture event if the companion 
started at a distance sufficiently nearby for the AGB star to en¬ 
gulf it. The amount of orbital energy released is 


GMi,,M 2 GM1M2 
lai 2a 


2 X 10^^ erg, 


( 4 ) 


where Mii = \.6Mq (assuming that the present circumstellar 
material of 1 M© comes from the AGB star), and assuming 
( 22 ,i = 4 AU (sufficiently near to be captured). 

There are considerable uncertainties in the orbital and energy 
estimates, and the common-envelope evolution sc enario is com¬ 
plex w ith an uncertain energy transfer efficiency (llvanova et al.l 
I 2 OI 3 I) . But there is good evidence that we see the effect of a 
capture event in which a fair fraction of the mass of the AGB 
star was released as the companion spiraled inwards. Most of it 


remains in the equatorial plane of the binary, but some of it has 
been accelerated up to high velocities in the form of a bipolar jet. 
The complex structure at the end of the hour-glass structure may 
be the imprint of what took place during the capture. However, 
taken at face value, the released orbital energy is not enough, and 
another mechanism may augment, or even dominate, the ejec¬ 
tion event, for insta nce a magnetic field o f the (rapidly) spinning 
[v sin / 50km s“^, IBakker et~^ (Il996ah ] post-AGB star. 

4. Conclusions 

We have presented CO isotopologue emission line data on the 
post-AGB object HD 101584 obtained with ALMA and pro¬ 
vided an initial analysis of the morphology, kinematics, masses, 
and energetics of the circumstellar medium. We conclude that 
the circumstellar medium of HD 101584 has been severely af¬ 
fected by the evolution of the binary system. A high-velocity jet 
(^ 150kms“^), directed at PA 90° and slightly tilted w.r.t. the 
line of sight, has excavated an hour-glass structure in the cir¬ 
cumstellar medium. A complex structure is seen at the extreme 
velocities of the hour-glass structure. The estimated ages of the 
complex structure and the jet are comparable, 500 yr. A com¬ 
mon origin of these phenomena such as the infall of the com¬ 
panion is therefore possible. However, taken at face value, the 
kinetic energy of the accelerated gas exceeds the released orbital 
energy during the capture event, and, further taking into account 
the expected efficiency of the process, the observed phenomenon 
does not match current common-envelope scenarios. This sug¬ 
gests that another process must augment, or even dominate, the 
ejection process. There appears to be substantial amount of ma¬ 
terial, maybe as much as a solar mass, in the central region, 
which is characterized by narrow spectral features that are only 
a few km s“^ wide. This may be the gas removed from the AGB 
star when the companion spiraled inwards. At low brightness 
levels there are indications of another bipolar structure. Its ve¬ 
locity gradient is opposite to that of the CO outflow and directed 
at PA -45°. 
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